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Abstract 

A truncation scheme for the Dyson-Schwinger equations of QCD in 
Landau gauge is presented which implements the Slavnov-Taylor iden- 
tities for the 3-point vertex functions. Neglecting contributions from 4- 
point correlations such as the 4-gluon vertex function and irreducible scat- 
tering kernels, a closed system of equations for the propagators is obtained. 
For the pure gauge theory without quarks this system of equations for the 
propagators of gluons and ghosts is solved in an approximation which al- 
lows for an analytic discussion of its solutions in the infrared: The gluon 
propagator is shown to vanish for small spacelike momenta whereas the 
ghost propagator is found to be infrared enhanced. The running coupling 
of the non-perturbative subtraction scheme approaches an infrared stable 
fixed point at a critical value of the coupling, a c ~ 9.5. The gluon propa- 
gator is shown to have no Lehmann representation. The gluon and ghost 
propagators obtained here compare favorably with recent lattice calcula- 
tions. Results for the quark propagator in the quenched approximation 
are presented. 
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1 Introduction 



Despite the remarkable success of perturbative QCD the description of hadronic 
states and processes based on the dynamics of confined quarks and gluons re- 
mains the outstanding challenge of strong interaction physics. Especially, one 
has to explain why only hadrons are produced from processes involving hadronic 
initial states, and that the only thresholds in hadronic amplitudes are due to 
the productions of other hadronic states. To this end one would like to under- 
stand how singularities appear in the Green's functions of composite hadron 
fields where, on the other hand, they have to disappear in colored correlations 
functions. 

To study these aspects of QCD amplitudes non-perturbative methods are 
required, and, since infrared divergences are anticipated, a formulation in the 
continuum is desirable. Both of these are provided by studies of truncated sys- 
tems of Dyson-Schwinger equations (DSEs), the equations of motion of QCD 
Green's functions. Typically, for their truncation, additional sources of informa- 
tion like the Slavnov-Taylor identities, entailed by gauge invariance, are used 
to express vertex functions in terms of the elementary two-point functions, i.e., 
the quark, ghost and gluon propagators. Those propagators can then be ob- 
tained as selfconsistent solutions to non-linear integral equations representing 
a closed set of truncated DSEs. Some systematic control over the truncating 
assumptions can be obtained by successively including higher n-point functions 
in selfconsistent calculations, and by assessing their influence on lower n-point 
functions in this way. Until recently all solutions to truncated DSEs of QCD in 
Landau gauge, even in absence of quarks, relyed on neglecting ghost contribu- 
tions completely §, |, |, f§. 

In addition to providing a better understanding of confinement based on 
studies of the behavior of QCD Green's functions in the infrared, DSEs have 
proven successful in developing a hadron phenomenology which interpolates 
smoothly between the infrared non-perturbative and the ultraviolet perturba- 
tive regime [|J, for recent reviews see, e.g., || In particular, a dynamical 
description of spontaneous breaking of chiral symmetry from studies of the DSE 
for the quark propagator is well established in a variety of models for the glu- 
onic interactions of quarks |H] . For a sufficiently large low-energy quark-quark 
interaction quark masses are generated dynamically in the quark DSE in some 
analogy to the gap equation in superconductivity. This in turn leads naturally 
to the Goldstone nature of the pion and explains the smallness of its mass as 
compared to all other hadrons. In this framework a description of the different 
types of mesons is obtained from Bethe-Salpeter equations for quark-antiquark 
bound states || . Recent progress towards a solution of a fully relativistic three- 
body equation extends this consistent framework to baryonic bound states, see 
e.g. flCfl and references therein. 

Here a simultaneous solution of a truncated set of DSEs for the propagators 
of gluons and ghosts in Landau gauge is presented [O, p"2| . An extension of this 
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self consistent framework to include quarks is subject to on-going research Jl3[ | . 
Preliminary results for the quark propagator in the quenched approximation 
have been obtained and will be shown. The behavior of the solutions in the 
infrared, implying the existence of a fixed point at a critical coupling a c « 9.5, 
is obtained analytically. The gluon propagator is shown to vanish for small 
spacelike momenta in the present truncation scheme. This behavior, though in 
contradiction with previous DSE studies @, |, |, [§, can be understood from 
the observation that, in our present calculation, the previously neglected ghost 
propagator assumes an infrared enhancement similar to what was then obtained 
for the gluon. In the meantime such a qualitative behavior of gluon and ghost 
propagators is supported by investigations of the coupled gluon ghost DSEs 
using bare vertices jl4|, |l{|. As expected, however, the details of the results 
depend on the approximations employed. 



2 The set of truncated gluon and ghost DSEs 

Besides all elementary 2-point functions, i.e., the quark, ghost and gluon prop- 
agators, the DSE for the gluon propagator also involves the 3- and 4-point 
vertex functions which obey their own DSEs. These equations involve succes- 
sively higher n-point functions. The gluon equation is truncated by neglecting 
all terms with 4-gluon vertices. These are the momentum independent tadpole 
term, an irrelevant constant which vanishes perturbatively in Landau gauge, 
and explicit 2-loop contributions to the gluon DSE. For all details regarding 
this truncation scheme we refer the reader to [12]. 

The ghost and gluon propagators are parameterized by their respective 
renormalization functions G and 

n m G ( fc2 ) n m fx KK\Z{&) 

D G\ k ) = p- ' D ^ k ) = [S^ - -p- J — jfea— • W 

In order to arrive at a closed set of equations for the functions G and Z , we 
use a form for the ghost-gluon vertex which is based on a construction from 
its Slavnov-Taylor identity (STI) which can be derived from the usual Becchi- 
Rouet-Stora invariance neglecting irreducible 4-ghost correlations in agreement 
with the present level of truncation Jl^] . This together with the crossing sym- 
metry of the ghost-gluon vertex fully determines its form at the present level 
of truncation: 

GM = iq »^) +W »\-GW)~ )' { ) 

With this result, we can construct the 3-gluon vertex according to procedures 
developed and used previously [|l6| , for details see |0| . 

We have solved the coupled system of integral equations of the present trun- 
cation scheme numerically using an angle approximation. The infrared behavior 
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of the propagators can, however, be deduced analytically. To this end we make 
the Ansatz that for x := k 2 — > the product Z(x)G(x) — > cx K with k ^ and 
some constant c. The special case k = leads to a logarithmic singularity for 
x — > which precludes the possibility of a selfconsistent solution. In order to 
obtain a positive definite function G(x) for positive cc from an equally positive 
Z(x), as x — > 0, we obtain the further restriction < k < 2. The ghost DSE 
then yields, 

CM - (>7? (i - i))"' c-'x- * * W ( 9 > 7? (I - 1)) cV" , 

where 7g = 9/(647r 2 ) is the leading order perturbative coefficient of the anoma- 
lous dimension of the ghost field. Using these relations in the gluon DSE, we find 
that the 3-gluon loop contributes terms ~ x K to the gluon equation for x — ► 
while the dominant (infrared singular) contribution arises from the ghost-loop, 

Requiring a unique behavior for Z(x) we obtain a quadratic equation for k with 
a unique solution for the exponent in < k < 2: 



^61-71897 ^0.92. (3) 
19 v ' 

The leading behavior of the gluon and ghost renormalization functions and thus 
of their propagators is entirely due to ghost contributions. The details of the 
approximations to the 3-gluon loop have no influence on the above consider- 
ations. Compared to the Mandelstam approximation, in which the 3-gluon 
loop alone determines the infrared behavior of the gluon propagator and the 
running coupling in Landau gauge (l], |[ ||, [|, this shows the importance of 
ghosts. The result presented here implies an infrared stable fixed point in the 
non-perturbative running coupling of our subtraction scheme, defined by 

a s (s) = fz(.s)G\.s) - ^-(--l) '^9.5, (4) 



4tt v ' x ' 9 V K 2 , 

for s -> 0. This is qualitatively different from the infrared singular coupling of 
the Mandelstam approximation B. 



3 Comparison to lattice results 

It is interesting to compare our solutions to recent lattice results available for 
the gluon propagator |l7j] and for the ghost propagator [M using lattice versions 
to implement the Landau gauge condition. We would like to refer the reader 
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to ref. [|19| where this has been done in some detail. It is very encouraging 
to observe that our solution fits the lattice data at low momenta rather well, 
especially for the ghost propagator. We therefore conclude that present lattice 
calculations confirm the existence of an infrared enhanced ghost propagator 
of the form Dq ~ l/(fc 2 ) 1+K with < K < 1. This is an interesting result 
for yet another reason: In the calculation of |Q the Landau gauge condition 
was supplemented by an algorithm to select gauge field configurations from the 
fundamental modular region which is to avoid Gribov copies. Thus, our results 
suggest that the existence of such copies of gauge configurations might have 
little effect on the solutions to Landau gauge DSEs. 

Here we want to add a remark concerning the comparison of the running cou- 
pling obtained in our calculation to lattice results. Recent lattice calculations of 
the running coupling are reported in Rcfs. [pp| based on the 3-gluon vertex, 
and Ref. ||| on the quark-gluon vertex. The non-perturbative definitions of 
these couplings are related but manifestly different from the one adopted here. 
One of the most recent results from the 3-gluon vertex is shown in the left graph 
of Fig. Q and compared to the three-loop expression which is for the momenta 
displayed almost identical to our expression (4) for the running coupling. This 
lattice result is obtained from an asymmetric momentum subtraction scheme. 
This corresponds to a definition of the running coupling g~l GVas which can ex- 
plicitly be related to the present one (g 2 (t,g) with t = ln////^ and g := g(p)), 

-2 U 2, - v . 2, v G 2 (.s) / P(g(t,g)) V 

g (t,g )3GVas = g (t,g ) km 2 1 _. , . (5) 

An inessential difference in these two definitions of the running coupling is the 
last factor in brackets in eq. (|^) which can be easily accounted for in compar- 
ing the different schemes. However, the crucial difference is the ratio of ghost 
renormalization functions G(s — > 0)/G(/i' 2 ). These considerations show that 
the asymmetric scheme can be extremely dangerous if infrared divergences oc- 
cur in vertex functions as our calculation indicates. Clearly, from the infrared 
enhanced ghost renormalization function this scale dependence could account 
for the infrared suppressed couplings which seem to be found in the asymmetric 
schemes. 

Similarly, the results from the quenched calculation of the quark-gluon ver- 
tex of Ref. [E2I which are compared in the right graph of Fig. [j] to our solution 
are obtained from an analogous asymmetric scheme. It is thus expected to have 
the same problems in taking the possible infrared divergences of the vertices 
into account which arise in both, the 3-gluon and the quark-gluon vertex, as a 
result of the infrared enhancement of the ghost propagator. 

Furthermore, definitions of the coupling which lead to extremas at finite val- 
ues of the scale correspond to double valued /3-functions with artificial zeros. If 
the maxima in the couplings of the asymmetric schemes at finite scales are no 
lattice artifacts, these results seem to imply that the asymmetric schemes are 
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P=6.2 ; Volumc=24 




1.0 2.0 3.0 4.0 n[GeV] 

Figure 1: Lattice results of the running coupling from the 3-gluon vertex (left, 
together with a 3-loop fit, Fig. 1 of Ref. |21|]), and from the quark-gluon vertex 
for (3 = 6.0 on a 16 3 x 48 lattice (right, c.f. Fig. 2 of Ref. Q). 



less suited for a non-perturbative extension of the renormalization group to all 
scales. Indeed, the results for the running coupling from the 3-gluon vertex ob- 
tained for the symmetric momentum subtraction scheme in Ref. pl| ] differ from 
those of the asymmetric scheme, in particular, in the infrared. These results 
would be better to compare to the DSE solution, however, they unfortunately 
seem to be much noisier thus far (see Ref. plf ). 

The ultimate lattice calculation to compare to the present DSE coupling 
would be obtained from a pure QCD calculation of the ghost-gluon vertex in 
Landau gauge with a symmetric momentum subtraction scheme. This is unfor- 
tunately not available yet. 



4 Quark Propagator 

We have solved the quark DSE in quenched approximation juj. In a first step we 
have specified the quark-gluon vertex from the corresponding Slavnov-Taylor 
identity. It contains explicitely a ghost renormalization function, 

r^p,q)=G(k 2 )r^(p,q) (6) 

where T^p is the Curtis-Pennington vertex (for its definition see e.g. |7j). It 
is obvious that this leads to an effective coupling very different from the one 
in Abelian approximation, especially in the infrared: This effective coupling 
vanishes in the infrared and is similar to the lattice result of Ref. p2|| shown in 
the right graph of Fig. [j] the main difference being that the maximum occurs 
at lower scale, /x ss 220MeV. This leads to a kernel in the quark DSE which is 
only very slightly infrared divergent. This allows, e.g., to use the Landshoff- 
Nachtmann model for the pomeron in our approach. With our solution we 
obtain as Pomeron intercept 2.7/GeV as compared to the value 2/GeV deduced 
from phenomenology, see e.g. pq]. 
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We have found dynamical chiral symmetry breaking in the quenched ap- 
proximation. Using a current mass, m(lGeV) = 6MeV we obtain a constituent 
mass of approximately 170 MeV. In the Pagels-Stokar approximation the cal- 
culated value for the pion decay constant is 50 MeV. These numbers are quite 
encouraging, especially for proceeding with the self-consistent inclusion of the 
quark DSE into the gluon-ghost system. 

Considering the quark loop in the gluon DSE one realizes that the quark loop 
will produce an infrared divergence which is, however, subleading as compared 
to the one generated by the ghost loop. In the latter there appear three ghost 
renormalization functions in the numerator and one in the denominator leading 
effectively to an infrared divergence of the order (fc 2 )~ 2K . In the quark loop term 
there is only one factor G and thus a divergence of type (k 2 )~ K . Due to this 
subleading divergence the infrared analysis has to be redone completely before 
one is able to draw conclusions whether or not and how quark confinement is 
implemented in our set of truncated DSEs. 

5 Summary 

In summary, we presented a solution to a truncated set of coupled Dyson- 
Schwinger equations for gluons and ghosts in Landau gauge. The infrared 
behavior of this solution, obtained analytically, represents a strongly infrared 
enhanced ghost propagator and an infrared vanishing gluon propagator. 

The Euclidean gluon correlation function presented here can be shown to 
violate reflection positivity jl2j , which is a necessary and sufficient condition for 
the existence of a Lehmann representation. We interpret this as representing 
confined gluons. In order to understand how these correlations can give rise to 
confinement of quarks, it will be necessary to redo the infrared analysis including 
self-consistently the quark propagator. Nevertheless, we found dynamical chiral 
symmetry breaking in the quenched approximation. 

The existence of an infrared fixed point for the coupling is in qualitative 
disagreement with previous studies of the gluon DSE neglecting ghost contri- 
butions in Landau gauge jl[ |5[ ||[ @]. On the other hand, our results for the 
propagators, in particular for the ghost, compare favorably with recent lattice 
calculations jl7|, [18| . This shows that ghosts are important, in particular, at 
low energy scales relevant to hadronic observables. 
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